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Abstract 

We review the current state of photon and dilepton measurements at RHIC, emphasizing that of the theoretical 
work seeking to interpret them. We highlight the progress made recently in the modelling of relativistic nuclear 
collisions, and explore the effect on electromagnetic observables. Some outstanding puzzles are presented. 
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1. Introduction 

The behaviour of strongly-interacting matter is known to be dictated by QCD, and colliding large nuclei at rel- 
ativistic energies has been convincingly shown to be an efficient means of studying QCD under extreme conditions. 
A vibrant experimental program is currently under way at RHIC (the Relativistic Heavy Ion Collider, at Brookhaven 
National Laboratory), and at the LHC (the Large Hadron Collider, in CERN). In the analyses of experiments per- 
formed at those facilities, measurements of the produced hadronic particles have revealed many fascinating features 
of the relativistic many-body systems. Whereas energetic jets contain information on the nature of the QCD plasma at 
early times Q]|2), me softer hadrons have shown an impressive final-state collectivity that has become a cornerstone 
of the qualitative success of hydrodynamics at RHIC and LHC energies Ql|4|. Electromagnetic observables com- 
plement well the hadronic measurements. Owing to the relative smallness of the fine-structure constant, a, real and 
virtual photons emerge unscathed from their production site: They probe the entire space-time trajectory of the nu- 
clear collision, from hot QCD plasma to the cooler hadronic phase. The price to pay for such a rich probe is a modest 
emission rate and a small signal-to-background ratio, in what concerns thermal photons, real and virtual. We review 
here some recent progress in the calculation of thermal photon yield and flow, and highlight some of the puzzles that 
have emerged in the analysis of electromagnetic radiation in relativistic heavy ion collisions. 



2. Electromagnetic emission rates 

Throughout the history of a single nuclear collision event, several sources of photons will come into play. The 
current formulation of relativistic hydrodynamics assumes that thermalization sets in after a proper time to. The 
value of this parameter is currently not calculable with certainty: it is left as a variable. Consequently, no reliable 
estimates exist for pre-equilibrium photons, other than those being generated in the very first instants of the collisions, 
calculable with perturbative QCD (pQCD), and measured in pp collisions. These "prompt photons" constitute an 
irreducible background to the photons being generated through some influence of a thermal medium. This statement, 
however, should be clarified: pQCD photons have a jet fragmentation component which may be sensitive to quenching 
effects from the medium [5 1 . The QCD jets interacting with the medium then convert into photons [6 7|. Atr > tq, the 
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thermal components of the medium may also interact to generate photons. There, the appropriate respective emission 
rates, R, for real and virtual photons can be expressed using finite-temperature field theory (8) as 

cP R g^ v 1 

w 5% = -W^T ImI1 " (w ' k) ' 
E+E -^k: = ^^[^- + ^--^-H^n>,k)-^ (i) 

In the above, lepton masses have been neglected, Im n R is the retarded photon self-energy, a> and k are components of 
the photon four-momentum k, E ± and p ± are components of the lepton/antilepton four-momenta, and ft = l/T. Those 
expressions receive correction of 0{a 2 ) in the electromagnetic interaction, but are exact in the strong interaction. 

The finite-temperature photon emission rates in the quark-gluon plasma (QGP) phase have been worked out to 
order a s , the leading order in the strong interaction, a little more than a decade ago ll9 fflD . One may question the 
robustness of those rates against NLO corrections, especially in light of closely related investigations of the behaviour 
of higher-order correlation functions lfT2H T3 l. This has been done recently, and first results have been reported [14|. 
It is shown there, that the phenomenological consequences of the NLO photon-emission rates are but a modest ~2Q% 
enhancement in the region of momentum |k|/T w 10. It is inappropriate here to discuss the technical intricacies of this 
calculation, but suffice it to say that it therefore appears that conclusions reached through previous analyses relying 
on LO QGP rates need not be amended. From the theoretical point of view, however, this work represents definite 
progress. 




Figure 1: The rate for the production of back-to-back lepton pairs, calculated in quenched lattice QCD at T ^ 1.45 T c , for two flavours of quarks. 
The top two curves of the legend represent two different ways to extract the lattice data. Also shown are calculation done with the Hard Thermal 
Loop resummation 1 15], and with keeping only the Born term. This plot is from Ref. [16]. 

In what concerns dileptons, another recent element of progress has been the nonperturbative evaluation of the 
vector current-current correlation function in lattice QCD at finite temperature. The connection with the dilepton 
emission rate is obtained by straightforward manipulations of Eq. ([T]l: 

d 4 R _ a 1 j n)ll 1 

where p n ' is the "normal" (8 | spectral representation of the current-current correlator, and M 2 = (p+ + p-) 2 . Such ex- 
trapolations from lattice calculations however require solving an inversion problem that is mathematically ill-defined: 
The lattice calculates Euclidian correlation functions, /(r), formulated in imaginary time, t. The connection with the 
physical Minkowski correlation function /(w), is 

/(t)) = dajf(a)) . ... — (3) 

Jo sinh (JSoj/2) 

Progress in inverting this (possibly noisy) integral has been obtained in recent years using statistical approaches like 
the Maximum Entropy Method IfTTl . Some some recent work has involved using a physically-motivated ansatz for the 
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lattice spectral function lfl"6l . The dilepton production rate obtained in this fashion is shown in Figure [T] The lattice- 
extracted dilepton production rate merges with the Born rate at cj/T ~ 4 and rises well over that, at low energies. 
Note that the low-M HTL rate grows differently as the invariant mass is lowered: A known behaviour iTTSl [T9l . 
More results of this nature, but for a range of temperatures closer to the transition temperature and at finite three- 
momentum, will eventually open the door to a precise non-perturbative characterization of the dilepton emission rates 
in the QGP domain. On the hadronic side, photon and dilepton rates have mostly relied on using effective hadronic 
Lagrangians, supported by the existing data on hadronic radiative decays, nuclear photoabsorption, and electron- 
positron annihilation [20-23 1. The different hadronic medium electromagnetic emissivities are summarized in the 
cited literature reviews: Current approaches do differ in detail, but agree on the importance of broadened spectral 
densities driven by in-medium interaction. 

Thus, there has been advances in the calculation of photon and dilepton emission rates. In addition to those, the 
field of relativistic heavy-on collisions has witnessed impressive leaps forward in the space-time modelling of the 
colliding system. As the electromagnetic rates need to be integrated to produce yields in order to be compared with 
experimental measurement, the space-time evolution needs to be considered on an equal footing. 



3. Space-time modelling 

As already mentioned, the success of relativistic hydrodynamics in the theoretical interpretation of relativistic 
nuclear collisions at RHIC and at the LHC is striking. In addition to momentum distributions, particle spectra may be 
characterized by their dependence on the azimuthal angle <p. This in turn will determine the momentum anisotropy, 
which can be quantified through a Fourier expansion [24]. The coefficients of this expansion, the flow coefficients, 
and the associated reaction plane are defined event-by-event: 

v„ = <cos (it (0 - ¥„))> , *F„ = - arctan . <Sm( "^ (4) 

n (cos (rup)) 

The recent progress in this area over the last few years has been such that the possibility of extracting transport 
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Figure 2: The photon Raa is shown, as a function of the photon momentum, measured in Au + Au collisions at RHIC. The curves shown are results 
from different theoretical approaches. Those are discussed in the text and in the references. The plot is from Ref. 1251 . 

coefficients from relativistic heavy ion data is a real possibility. More specifically, a small but finite value of the 
shear viscosity to entropy density density ratio (rj/s) appears consistent with a global fit to RHIC and LHC data sets, 
including measurements of higher flow harmonics [26|. 

Given the electromagnetic sources and emissivities discussed previously, what are the yields predicted by the 
modern hydrodynamical approaches? In what concerns real photons, a recent update [27] of data analysis by the 
PHENIX collaboration is shown in Figure [2] What is plotted there is photon Raa at RHIC: the ratio of the measured 
yields in Au + Au collisions, divided by the production cross section in p + p collisions multiplied by the average 
nuclear thickness function for the appropriate centrality. The theoretical curves of Figure prelate to calculations that 
included effects that stem from cold nuclear matter (isospin and (anti)shadowing), effects of photon generation by 
QGP and suppression by jet energy loss, and coherent effects such as LPM |25|. As shown in Figure [2] the data 
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currently can't distinguish between an approach with no QGP [28 1 and one which relies on 2+1D hydrodynamics 
where the photon-enhancing channels are almost canceled by the energy-loss effects [5 29]. Finally, an analysis of 
the low-momentum photon data, extracted from continuing dilepton measurements to the light-cone, has produced 
an effective temperature of the photon signal that exceeded pQCD expectations |30|. For central collisions: T = 
221 + 19 stat + 19 syst MeV. Even if the spread in initial temperatures in the theoretical models that seek to explain the 
same data is considerable, 300 < 7init(MeV) < 600, a consistent feature is that those values all exceed the typical 
transition temperature predicted by lattice QCD, as does the empirical value extracted from data. A critical review of 
the theoretical calculations that led to this wide palette of initial temperatures is due, and should reduce this scatter 
considerably. 

Dilepton spectra have been available at RHIC for some years. The first measurements there were done by the 
PHENIX collaboration. A striking feature of these is a large enhancement over background in the low invariant mass 
region: < M(GeV/c 2 ) < 0.7, in Au + Au collisions [ 30 1 . The enhancement factor over the hadronic cocktail 
expected from p + p measurements has been measured to be 4.7 + 0.4 stat + 1.5 syst + 0.9 model in minimum bias events, 
and this enhancement is concentrated at low pripr <1 GeV). Since these pioneering data were recorded, it is fair to 
write that they have constituted a challenge for theory. This mismatch between theory and experiment persists to this 
day - models underpredict the data typically by a factor ~ 3 at M ~ 400 MeV - and currently represents a puzzle. 




Figure 3: The yield of low mass dilepton as measured at RHIC by STAR (left panel) and by PHENIX (right panel). Also shown is the data to 
cocktail ratio, for the two respective experiments. The left panel is from Ref. 1311 . and the right panel from Ref. 1251 . 

Another element of experimental progress in the field has been the entrance of the STAR collaboration on the 
dilepton scene, owing to the recent installation of a Time-of-Flight detector. Both the TOF and the Time Projection 
Chamber (TPC) enable electron identification from low to intermediate pj. The recent STAR dilepton measurements 
are shown in Figure [3] This collaboration also reports an enhancement (over the expected hadronic cocktail) in the 
low dilepton invariant mass region. However, this enhancement is smaller than the corresponding PHENIX figure: the 
STAR factor is 1.53 + 0.07 stat + 0.41 syst . Results are still preliminary, but it appears unlikely that this difference in the 
enhancement factor observed by the two different collaborations is due to acceptance effects l32l : it is thus currently 
also a puzzle. Much is expected in the low invariant mass region from analyses already made with PHENIX's Hadron 
Blind Detector, designed to reduce background by suppressing Dalitz pairs and photon conversions ||33| . 

3.1. Viscous effects 

The last few years have seen the advent of fully three-dimensional hydrodynamics simulations, with viscous 
corrections up to second order in flow velocity gradients: see the previously cited reviews ((3] |4)), an d references 
therein. We will only discuss here the effect of shear viscosity; How will it affect the generation of electromagnetic 
radiation? The introduction of a finite value of the shear viscosity coefficient will modify the microscopic distribution 
functions. In a multispecies ensemble, a popular ansatz for modifying the equilibrium distribution function of species 
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i is 

foi -> h + Sfi , Sfi = / 0i (l + fodp a pPx af } — * „ 2 , (5) 

where n a p is the shear correction to the stress-energy tensor, e is the energy density, and P is the pressure. This ansatz 
is however not unique l34l . In the QGP phase, the complete leading -order rates in a s have yet to be re-evaluated, 
but the viscous corrections to the Compton and annihilation contributions have been done [35-38]. In the hadronic 
(confined) sector, the evaluation of viscous corrections to photon emission rates requires a recalculation of the many 
different channels that went into the compilation of Ref. [22], for example, using Eq. Q. The viscous corrections 
to the photon yields thus manifest themselves through a modified hydrodynamic evolution, and also by altering the 
microscopic rate equations themselves [38). The respective and combined effects are shown in Figure[4j for conditions 
corresponding to Au + Au collisions at -\fs = 200 GeV, and a - 10% centrality class. For the partonic phase photons, 




Figure 4: The photon yields, after integrating the rates through the hydrodynamic evolution. The full curves are results obtained with inviscid 
hydrodynamics, the dotted lines represent yields obtained with viscous hydro, but using uncorrected rates, and the dashed-dotted curves are yields 
with the viscous hydro and with the appropriately modified emission rates. A value of r//s = 1/4jt has been used. Left panel: Photons produced 
from the parton phase only. Right panel: Photons from the hadronic gas only. This plot is from Ref. (38J. 

using the viscous hydro (without changing the rates) reduces the yield at high pj\ the viscous hydro requires a lower 
initial temperature (than ideal hydro) to be compatible with hadronic data, as entropy is generated throughout the 
evolution. For the hadronic gas photons, the main effect of viscosity is to produce a flatter temperature-against-time 
evolution [38], and the higher late-stage temperatures will then produce more photons. From Eq. Q, the viscous 
corrections to the rates are proportional to momentum and to the shear pressure tensor. The tensor components are 
maximal at small times and high temperatures, driven by the large initial gradients, and are much smaller at later 
times [ 38 1 . This supports the net viscous yields in the hadronic gas phase not being noticeably affected by the viscous 
corrections to the microscopic processes. Conversely, the partonic photons receive a correction which grows with 
momentum, making the spectra harder [36, [37] [39). However, this hardening of the spectra is modest, building up to 
a factor of m 2 at pj ~ 4 GeV. The pQCD photons will take over in this momentum region, making the extraction of 
a viscosity from photon spectra a challenging task. 

3.2. Fluctuating initial conditions 

Another recent improvement in the understanding of nuclear dynamics, is that of the role payed by initial con- 
ditions and of their effect on finite, odd flow harmonics BUI for hadrons. Including Fluctuating Initial Conditions is 
therefore now an unavoidable feature in the modelling of relativistic nuclear collisions. How will FICs affect the pro- 
duction of electromagnetic radiation? One must first define the nature of the FICs: Here we report on results obtained 
with a Monte Carlo Glauber approach, previously used for making a successful prediction of the hadronic V3 at RHIC 
ETll . The lumpy initial states germane to FICs lead to harder hadron spectra, owing to a combination of hot spots 
and of larger blue-shifts from initial gradients J42). The same arguments will apply to photons, see also Ref. [43 1. 
On Figure [5] we show the effects of FICs, combined with those of a finite shear viscosity coefficient. An immediate 
conclusion from that figure: the modification to the spectra from either a finite 77 or from FICs are of similar magni- 
tude. As hydrodynamics is inherently non-linear, it is thus important to have both ingredients in a single calculation. 
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Figure 5 : Left panel: Photons from the QGP, using averaged initial conditions for the relativistic hydrodynamics evolution (ideal hydro: solid curve; 
viscous hydro: dash-dotted curve), and using FICs (with viscosity): double-dotted curve. Right panel: Same as left panel, but but for hadronic gas 
photons. This plot is from Ref. [38]. 



Again, the effects are largest in regions where pQCD photons will set in. From a theoretical point of view, to make the 
electromagnetic emissivity consistent with advances in hadronic modelling should however be interpreted as a clear 
element of progress. 

4. Photon elliptic flow 

As the measurement of higher flow harmonics have contributed significantly to the success of relativistic hydro- 
dynamics by quantifying the flow anisotropy, it is natural to inquire about the power of similar observables in the 
electromagnetic sector. As real and virtual photons decouple from the matter once they're produced, and because the 
photon production rates are mostly f-channel dominated (depending on the transverse momentum of the produced 
photon), their elliptic flow should track the matter momentum anisotropy from early (high-p^ photons) to late (lower- 
Pt photons) times fl4l . Thus, photon elliptic flow is therefore a discriminating observable, complementary to parallel 
hadronic measurements. A similar story is predicted to hold true for dileptons [45 1. It is also important to investigate 
the consequence of out-of-equilibrium effects, as quantified by a finite coefficient of shear viscosity, and of FICs on 
photon elliptic flow. These effects on v y 2 have been studied individually [36, 37], and together 1581 . They are shown 
in Figure [6] As for hadrons ATI , viscous effects reduce the photon momentum anisotropies: V2 is made smaller by a 
finite r\. In this centrality class, the FICs make the elliptic flow larger but they also could lead to a reduced V2 in more 
peripheral events ATI . As in the case of spectra, it is therefore important to treat viscous corrections and fluctuating 
initial states simultaneously, as they compete in magnitude. 

Preliminary photon v% data have been available for some years, but were limited to a small range in pj and to low 
precision [46, 47 1. This situation has been firmed up recently, and new results from the PHENIX collaboration have 
been released [48 1 . At higher pj, the direct photon measurements are consistent with a vanishing elliptic flow, but 
the experimental uncertainties are currently too large to rule out the small positive contribution from jet fragmentation 
and/or the small negative contribution from jet-photon conversions 11231 . At low p T however, the photon elliptic flow is 
positive, as is that of the strongly interacting source, but is surprisingly as large as that of the 7r° flow and increases with 
centrality. Current models based on relativistic hydrodynamics - which are successful in interpreting hadronic flow - 
predict much smaller values for photon elliptic flow, as is also clear from Figure [6] see also [49 1 . This experimental 
result has since been confirmed by an independent measurement of the conversion electrons [32|. 

One may start from the vZ measurements, and investigate the demands they put on the collisions dynamics. This 
aspect was pursued in Ref. Il50l . There, an elliptic fireball expansion constrained by, and consistent with, bulk hadronic 
data was used. This empirical fireball approach is also consistent with photon measurements at CERN energies 
performed by the WA98 collaboration, and with dilepton measurements performed by the NA60 collaboration. A 
configuration which appeared consistent with existing data yields the photon spectrum and elliptic flow shown in 
Figure [7] A key ingredient in this work is the need for a large acceleration of the fireball, which blue-shifts the photon 
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Figure 6: Left panel: Thermal photon vi from the QGP only. Right panel: Net photon elliptic flow. The solid curve is the result of using inviscid 
hydrodynamics with averaged initial conditions. The dot-dashed curve shows the effect of a finite shear viscosity coefficient, and the double-dotted 
curve shows the combined effects of r\ and FICs. This plot is from Ref. [38). 



spectra to the required hardness. A consequence is then the reduced contribution of the QGP plasma radiation, and 
the domination of the Hadronic Gas phase, in the HG phase, flow has time to build up and the resulting photon 
elliptic flow is now large enough to clip the bottom of the experimental error bars. Note however that neither viscous 
corrections nor fluctuating initial conditions are explicitly considered here, other than through the effective dynamics. 
Further work is needed to reconcile these results with those obtained through detailed dynamical space-time models. 




p T (GeV) p T (GeV) 



Figure 7: Left panel: Thermal photon yields shown together with PHENIX data. Right panel: the related elliptic flow coefficient. The photons for 
the partonic phases (QGP) are show, and so are those from the hadronic gas (HG). The other curves are explained in the cited reference. This plot 
is from Ref. (50). 



5. Conclusions 

Measurements of electromagnetic radiation from colliding relativistic nuclei have reached an impressive degree 
of sophistication, even if puzzles remain in what concerns the consistency of different measurements, as well as in 
their theoretical interpretation. The theory effort has seen considerable consolidation in the modelling of hadronic 
observables and in the calculation of electromagnetic radiation. More work needs to be done, and this will continue. 
On the experimental side, next-generation results from RHIC as well as new results on photons and dileptons from 
the LHC will continue to drive the field forward. 
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